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Abstract 

At sufficiently high temperature and density, quantum chromodynamics (QCD) is expected 
to undergo a phase transition from the confined phase to the quark-gluon plasma phase. 
In the Lagrangian lattice formulation the Monte Carlo method works well for QCD at 
finite temperature, however, it breaks down at finite chemical potential. We develop a 
Hamiltonian approach to lattice QCD at finite chemical potential and solve it in the case 
of free quarks and in the strong coupling limit. At zero temperature, we calculate the 
vacuum energy, chiral condensate, quark number density and its susceptibility, as well 
as mass of the pseudoscalar, vector mesons and nucleon. We find that the chiral phase 
transition is of first order, and the critical chemical potential is fic = m iy*„ (dynamical 
quark mass at /i = 0). This is consistent with hq w Afjy /3 (where is the nucleon 
mass at \i = 0). 
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1 Introduction 



1.1 Motivation 

According to the big bang model in cosmology, the early universe underwent a series of 
drastic changes. For some time it was a hot and dense quark-gluon plasma (QGP), where 
quarks and gluons were deconfined. Today it is in a low temperature and low density 
hadronic phase, where quarks are confined. The ultimate goal of machines such as the 
Relativistic Heavy Ion Collider (RHIC) at BNL and the Large Hadron Collider (LHC) at 
CERN is to create the QGP phase. The QGP may also exist in the core of very dense 
stars such as neutron stars. Quantum chromodynamics (QCD) is the fundamental theory 
of quarks and gluons. A precise determination of the QCD phase structure at finite tem- 
perature T and chemical potential /x will provide valuable information in the experimental 
search for the QGP. Lattice gauge theory (LGT) proposed by Wilson in 1974, is a very 
reliable technique for the investigation of phase transitions. There are no free parameters 
in LGT when the continuum limit is taken, in contrast with other nonperturbative tech- 
niques. Although the standard lattice Lagrangian Monte Carlo method works very well 
for QCD at finite temperature, it unfortunately breaks down at finite chemical potential 
(due to the so-called complex action problem). This is briefly summarized in Sect. |1.2j . 
On the other hand, lattice QCD at finite chemical potential formulated in the Hamiltonian 
approach does not encounter a complex action problem. In Sect. 0, we develop a Hamil- 
tonian approach to lattice QCD at finite chemical potential \x. We solve this in the case of 
free quarks and in the strong coupling limit. 

1.2 Present status 

LGT is an approach to QCD from first principles. However, it is not free of problems: 
(a) First, there are lattice artifacts: A finite volume and a finite lattice spacing introduce 
errors, (b) There is a no-go theorem for chiral fermions: There is species doubling of any 
local fermionic theory with continuous symmetries. For naive fermions, chiral symmetry is 
preserved, but the species are doubled and the chiral anomaly is wrong. Kogut-Susskind 
fermions preserve the continuous U(l) chiral symmetry, but break explicitly flavor sym- 
metry. For Wilson fermions, the flavor symmetry exists, but chiral symmetry is explicitly 
broken. Kogut-Susskind fermions and Wilson fermions have been extensively used in nu- 
merical simulations. Recently, there has been evidence showing that those two approaches 
may give the topological charge or anomaly incorrectly ]I[ on a finite lattice. Therefore, 
it is far from clear whether correct results in the continuum can be obtained using those 
fermion formulations. Kaplan's domain wall fermions || and Neuberger's overlap fermion 
formulation || have attracted much attention, because they give the correct chiral modes, 
they also produce the correct anomaly and topological charge. For domain wall fermions 
there is an extra dimension and the lattice size in this dimension has to be very large. Thus 
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algorithms suitable for those new fermion approaches need to be developed. In this paper, 
we do not address those problems. 



Here we would like to investigate lattice QCD at finite chemical potential. In the contin- 
uum, the grand canonical partition function of QCD at finite temperature T and chemical 
potential /x is given by 

Z = Tr e -/3(H-„N) ; p = Qt B T)-\ (1) 

where ks is the Boltzmann constant, H is the Hamiltonian, and N is particle number 
operator 

N = [ d 3 x ^(x)ip(x). (2) 



The energy density of the system with free quarks is given by || 

e=--TrHe-^- N ) = -i^U. (3) 
V Z V d/3 v ; 

Going over to T — > 0, the energy density (where the contribution of /x = is subtracted) 
becomes 

4_7r r°° I 3 \ 3 

e sub = J ^ d 3 pQ U- , + m 2 I J] v) + m 2 . (4) 

Here is the step function. In the chiral limit m — > one obtains 

1 r°° 

6sub = 4^ J_ 00 d3 P Q( J i -\P^\^ 

- ?jTw*^ = ^. (8) 

In the Hamiltonian formulation of LGT, Eq.(^) in well defined. For Wilson fermions or 
Kogut-Susskind fermions, the relation Eq.(||) is satisfied (see below). However, if one 
constructs the fermionic lattice Lagrangian via Legendre transformation of the Hamiltonian, 
one cannot reproduce the continuum relation Eq.(|5|). Let us take the naive fermions as an 
example. The action obtained via Legendre transformation of H reads 

a 3 ±4 

S f = a A J2mij;(x)^(x) + — Y Y if)(x)^ k ^{x + k) + aVX)^ ( x )i>(x), (6) 

X ^ x k=±l x 

where = — 7^. This action gives the following result for the subtracted energy density 

-Q~ 4 r , 4 Dj=i sin2 Pj + ( mQ ) 2 r _ nl m 

esub 4tt 4 J-n P (sin p 4 -^ a ) 2 + £ 3 =1 sin 2 Pj + (ma) 2 [f * J ' 1 1 
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Taking the limit m — > and the continuum limit a — > 0, e su b oc (fi/a) 2 , i.e. becoming 
quadratically divergent, and therefore it is inconsistent with the continuum result of Eq.(H). 
This problem is not due to the species doubling of naive fermions, because the case of 
Kogut-Susskind fermions or Wilson fermions is similar. 

Hasenfratz and Karsch || proposed the following solution: If (sinp4 — ifi) 2 is replaced 
by sin 2 (p 4 — ifi) the continuum result Eq.([5[) is reproduced, except for a factor of 16. 
Correspondingly in the action, the chemical potential is introduced in the following way 

a 3 3 r- » l 

x ^ x j=l 

3 

+ —Y, [e^^h^ix + 4) - e-^ix + l)^(x)] . (8) 



X 



The chemical potential can be introduced analogously for KS as well as for Wilson fermions. 
Such treatment of the chemical potential is numerically feasible in the quenched approxi- 
mation (where the fermionic determinant det A is constraint to be 1, and quark loops are 
suppressed). However, there is evidence || that the quenched approximation produces an 
unphysical onset of the critical chemical potential at the value fic — M w (m ^ 0)/2, being 
in conflict with other theoretical predictions fic ~ M$ /3 (M$ is the nucleon mass at 
fi = and M 7T (m ^ 0) is the pion mass at finite bare quark mass m. A finite bare quark 
mass must be introduced in most of the numerical simulations). The unphysical onset of 
/ic is considered as a defect of the quenched approximation. 

For full QCD, the fermionic degrees of freedom have to be integrated out. In the measure 
occurs the fermionic determinant det A. For finite chemical potential det A becomes com- 
plex (complex action problem), which renders numerical simulations extremely difficult. 
Much effort has been made to solve the notorious complex action problem: 

(1) The Glasgow group has suggested to treat det A as observable This method requires 
a very large number of configurations, in particular for \i « fic- Even on a very small 
lattice V = 4 4 , the computational costs exceed the current computer capacity ||. 

(2) In the imaginary chemical potential method det A becomes real, which works well 
for numerical simulations at high temperature and low density. But it might not work at 
low temperature and high density. 

(3) It has been proposed to utilize a special symmetry Rl0| . This is the only successful 
method in Lagrangian lattice QCD, but it works only for the SU(2) gauge group. 

(4) Recently, a new approach has been proposed in pl]| , using quantum spin variables. It 
remains to be seen whether this can be applied to QCD. 
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2 Hamiltonian Approach 

2.1 Free fermions at zero chemical potential 

The lattice Hamiltonian describing noninteracting Wilson fermions in d + 1 dimensions at 
fi = reads 

_ j _ ~ f r _ _ ~ - 

H = mip(x)i()(x)-\ ip(x)'jk4 , {x+k)-\ \ip(x)%lj(x) — ip(x)ip(x + k) . (9) 



x,k=±j 



x,k=±j 



We want to diagonalize H so that the fermionic field ip can be expressed in terms of up 
and down 2-spinors £ and , 

5 • 



(10) 



We define the bare vacuum state |0 > as 

^|0> = 77|0> = 0. 



Since the up and down components are coupled via the 7^ matrices, the bare vacuum is 
not an eigenstate of H. Let \Q > denote the physical vacuum state, and Eq the vacuum 
energy. One can use a unitary transformation to decouple the up and down components 
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H' = exp(— iS) H exp(iS). 



(12) 



Such a transformation is similar to the Foldy-Wouthuysen transformation [|T^]. Then the 
physical vacuum state of H can be expressed as 

\Q) =exp(iS)|0). (13) 

The operator S can be computed explicitly. For Wilson (r 7^ 0) or naive (r = 0) fermions 



it reads 12 



S 



S ,, 



A, 



A p j=i 



a 

1/2 



E 

3=1 



(14) 
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and p is the momentum. The transformed Hamiltonian becomes 



2r 



£sin 2 (pj-a/2) 



cos 2# p + A p sin 2# p 



+ 



cos 261 



2r 



m 



+ -$>n 2 (p i a/2) 



sin 29 v 



ifjplfjp 

sinpjd 



Xl^Ti— /<• 



(15) 



The vacuum energy is given by 

E n = = (0\H'\0) 

= -2N c N f Y: 



2r 



m 



+ -E sin2 (pW 2 ) 



cos 29 p + A p sin 29 p 



(16) 



where iV c and Nj, respectively, are the number of colors and number of flavors. The vacuum 
energy Eq is minimized under variation of the parameters 9 P if 



tan 26L 



A, 



(17) 



This condition also leads to the cancellation of the second term in Eq. (|T|) coupling the 
up and down components such that 

H'\0) = Y,A' p ^p\0) = E n \0), (18) 



where we denote 



A' 



The vacuum energy becomes 



2r d 

m H ^2 sin 2 (Pj a /2) 

a 3=1 



E n = -2N c N f J2A' p . 



+ Al 



1/2 



(19) 



(20) 



It can be easily seen that \fl) is the eigenstate of H and Eq is its eigenvalue. For Wilson 
fermions, in the continuum limit a — > 0, for any finite momentum p, we have 



A' 



m 2 + p 2 , 



(21) 



giving the correct dispersion relation. 
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2.2 Free fermions at nonzero chemical potential 

We follow the same steps as in the case \x = 0. According to Eq. ([!]), the role of the 
Hamiltonian is now played by 

H^ = H- fiN, (22) 
where H is given by Eq.(|]) and N is given by Eq.(^). Let us define the state \n p ,n p ) by 

£p|0p;^p) 0) £p|0p; Tip) |lp3^p)> £p|lpj^p) |0p;^p); £p|lp;^p) 0, 

r) p \n p} P ) = 0, r?J|n p , P ) = \n p , l p ), r) p \n p , l p ) = \n p , p ), ?7j|n p , l p ) = 0. (23) 



The numbers n p and n p take the values or 1 due to the Pauli principle. By definition, 
the up and down components of the fermion field are decoupled. Obviously, this is not an 
eigenstate of due to the nondiagonal form of H. However, they are eigenstates of H' 
which are related to by a unitary transformation 



H'^ = exp(-iS) exp(iS) = H' - fxN. 
For the vacuum eigenstate of we make an ansatz of the following form 

= exp(iS)Y,fn p ,n p \n p ,n p ). 



(24) 



(25) 



S is given by Eq. (|T4"|) and the parameter 9 P is given by Eq . ([17|) . Both S and 9 P do not 
depend on \i because the quark number operator N commutes with S. H' is given by Eq. 
(|I5|) . The vacuum energy thus obeys 



E n = WH^\n) = J2fn p ,,n p Jn p n p (n p m p ,\H^\np,np) 

p',p 

= J2 C n P ,n P (n p ,n p \H' - fiN\n p ,n p ), 
p 

where we have introduced the notation C npfip = fn p ,n p - From Eq.(p3[) follows 

En = Yl C n P ,n p (A p (n p , n p \ip p ijj p \n p , n p ) - fi{n p , n p \i/)lif) p \n p , n p ) 



(26) 



2N c N f C np ,n p [(4 -fi)np+ (A' p + ^n p -A' p - fi 



(27) 



We have not yet specified the function C npt n p - For this purpose we use the condition of 
stability of the vacuum. Because \x > 0, the vacuum energy increases with n p . This means 
the vacuum is unstable unless n p = 0. This simplifies Eq. (|27D to 

(28) 



E n = 2N c N f £ C np [(A' p -^n p -A' p - \x 
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where we use the abbreviation C np = C np $. From the normalization condition Cq p +C\ p = 1, 
we obtain 

E n = 2N c N f £ [C lp (A' p - p) - A' p - p] . (29) 
p 

Cip depends on the value of p and its dependence can be seen by inspection of the derivative 

dEn 



2N c N f (A' p - p) . (30) 



For p > A' p , the right-hand side is negative. Maximizing C\ p means minimizing the vacuum 
energy. Therefore, C\ p = 1. On the other hand, for p < A' , the right-hand side is positive 
and for any C\ the vacuum is unstable. Therefore, C\ = 0. We can summarize these 
properties by writing 

C lp = 6 (p - A' p ) . (31) 

Thus the vacuum energy becomes 

En = 2N c N f Y / (C 1 A' p -A' p ). (32) 

p 

The subtracted energy density reads 

En-E n \v=o_ 2 ^ 2 3 



Here N s is the number of spatial lattice sites. In case of Wilson fermions, for m = and 
in the continuum a = 0, for any finite momentum p, one has A' = \p\. In 3+1 dimensions, 
at the corners of the Brillouin zone pja = (7r, 0,0), (0,7r, 0), (0,0, w), (tt, 7r,0), (0,7r, 7r), 
(7r, 0, 7r), (7r, tt, tt), one has 6 ^ — = 0. Therefore, in the continuum we find 

87T ft* u 4 



Thus we have proven that we can reproduce in the Hamiltonian formulation the continuum 
result of the vacuum energy density, Eq. (|5|). For naive fermions, in the continuum limit 
a = 0, there will be an extra factor of 2 d . 

2.3 Strong coupling QCD at nonzero chemical potential 
2.3.1 Structure of the Hamiltonian 

As is well known, lattice QCD at fi — confines quarks and spontaneously breaks chiral 
symmetry. For a sufficiently large chemical potential, this picture may change. At lattice 



spacing a ^ 0, as discussed in Sect. |1.2| , none of the standard approaches to lattice fermions 
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T A 


1 




74 


75 


^7475 


i747j 




%e i hi-fi <f1 hi h 


La 


1 


-1 + 25 k ,j 


-1 


-1 


1 


1 - 25 k j 


-1 + 25 k j 


1 - 2<5 fc j 



Table 1: T matrices and coefficients. 



is satisfactory. Here we set out to investigate finite density QCD in the strong coupling 
regime 1/g 2 << 1, using the Hamiltonian formulation. One of the goals is to get a better 
understanding of the mechanism of chiral phase transition. According to Ref. [0, H' in 
Eq. (p4|) now is replaced by 

(20 o )dl 



H' 



i-(26 yd 



+ 



a 



+ 



m 

9 C %aN^ ^ ^ ^i,/i( a: )7^ C 2,/i(a ; + ^)^ 2 ,/ 2 (a ; + ^)7fcV'c li / 2 (a;), (35) 



x k=±j 

where d = 3 denotes the spatial dimension, ci, c 2 are color indices and /i, / 2 are fla- 
vor indices (summation over repeated indices is understood), 9 = l/(4ma + g 2 C]y), and 
Cn = (N 2 — 1)/(2N C ). The four-fermion interaction is induced by gauge interactions with 
fermions. A very similar Hamiltonian has been derived in Ref. |14j] using strong coupling 
and large N c expansion. After a Fierz transformation, H' becomes [|12| 



H' 



m 



l-(29 ) 2 d 



+ 



(20 o ) d 



X 



g 2 C N d(26 

4a 



g 2 C N (2gg) 
32aN r 



L A ^ fl (x)T A ip f2 (x)ipj 2 (x + k)T A ip fl (x + k). 

x k=±j 



(36) 



The matrices T A and La are given in Tab.[l]. 
Let us define the following operators (12], [L5| , 

n(x)/ l/2 = 7r - r =^J l (x) (1 - 74) 75^/ 2 (^) 



2^1-v 

— A/ (y 



where 



1 



NfN s ~ 



^ ] C np ,n p (jlpi Tip I V'p V'p I ^p 5 ^p) 



2^ 

iV, 



^ Cn P ,n p ( n p + 



n, 



(37) 



(38) 



Using mean field approximation, one can show that fl5 



V 3 {x),V^x') =5 XjX , 



(39) 
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Thus the operators II and Vj, defined in Eq. (pTj), behave like pseudoscalar and vector op- 
erators. In Ref. |15| it has been shown that the operator ipif) satisfies the same commutation 
relations as v + 2ITTI + 2J2j VjVj. Therefore, H' in Eq. (0) can be written in terms of 
pseudo-scalar and vector particle operators in the following way 

+ V £ (W) 11 ^ + k ) + £ ^(^(z + *)(!- 2M + fc-c. ] , 



(40) 



where 



E. 



(o) 



«2 



l-(2^o) 2 ^ 



+ 



(20 o )d' 



m 



4a 



l-(29 ) 2 d 



y v v g 2 C N (29 yd , t _ 



(2# )d 



16aiV c 
g 2 C N d{29 ) 
AaN r 



-v. 



g 2 C N (29 ) 2 _ 
8aN c V ' 

£ C np ,n p (n p , n p \il)^ipp\rip, n p ) = E Cn v ,n v (n p — n p + 1) , 



V E ,% K > n p I $h jl>h ,P^h A , P \ n Pi n v) 



N f N, 
(2Nc 



f-L'S p 
2 



Cn p ,n p (n>p — n p + 1) , 

p 

tP \n p , Tip) 



NfN s — 



(41) 



In Eq . fl40j) , we have ignored the nonmeson terms which give no contribution to the energy. 
Making a Fourier transformation, one obtains 



h> = 4 0) + G 1 E(|n t b)n(p) + E^ t (p)^(p)j 

+ G 2 £ (n t (p)nt(-p) + /i.e.) £ cospj-a 
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+ G 2 E ( V j(P) V H~P) + E cosp^a - 2 cospja 



This can be diagonalized by a Bogoliubov transformation JT2] 



P(p) = cosh u(p)a(p) + sinhM(p)a^(— p), 
Vj(p) = cosh Vj(p)b(p) + sinh Vj(p)tf (— p), 



where 



tanh 2u(p) 
tanh 2vj(p) 



-2G 2 

gT 

-2G 2 
Gi 



E C0S Pi< 



cosjya — 2 cospj( 



(42) 



(43) 



(44) 



This condition also minimizes the vacuum energy. The Bogoliubov transformed Hamilto- 
nian eventually becomes 



z p 



1 - tanh 2 2u(p) - 1 



1 - tanh 2 2^(p) - 1 



+ G x E V 1 - tanh2 2M(p)a f (p)a(p) + d E V 1 " tanh2 2 ^(p) 6 J(p) 6 j(p)- ( 45 ) 



2.3.2 Vacuum energy 

The vacuum energy is given by 

E n = (n\H fl \n) = E n 0) ~N f N s ^ 



GiE 



1 - tanh 2 2u(p) - 1 



N 2 f 



+ ~y g ^ E v 1 - tanh2 2v M - 1 



■ (46) 



From Eqs. (|4i"D -(fi5D, we get for m = 



2N f N c 



E Gn p ,« p - p) n p + (m£ + p + 2m { ^ n n p ) n, f 

Tip,Tlp 



+ f g,e 

z p 



1 - tanh 2 2«(p) - 1 



N 



+ ^ G iE Vl-tanh^i;^)-! 

pj 



(47) 
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Here m^„ = d/(ag 2 CN) is the dynamical quark mass at n — 0. It is obvious that n p = 0, 
otherwise, the vacuum is unstable. Using the notation and normalization condition for the 



coefficient as in Sect. |2.2j , we obtain 



2N c N f 



+ 



iV 



1 - tanh 2 2u(p) - 1 



1 - tanh 2 2^- (p) - 1 



(4? 



Again, using the same argument as in Sect. p72| , the coefficient C\ must be 

C lp = 8 (ai - rog.) . (49) 

Substituting into Eq.([P|) yields 

go 
2JV«JV)iV. 



- /*) (/* ~ m dyn) - m dyn ~ H 



2iV. 



1 - tanh 2 2«(p) - 1 



2N. 



pj 



1 - tanh 2 2vj(p) - I5Q) 



2.3.3 Chiral condensate and critical // 

According to the Feynman-Hellmann theorem, the chiral condensate is related to the ground 
state energy by 



1 t <9E Q (m + 0) /T/xWh Q / (0) 



(51) 



where (ipip)^ is the chiral condensate at fj, — 

4ci \ 



(#) (0) -2iV c 1 



^C 2 



1-^-^7, 

iV c 1 iV c 



(52) 



and for d = 3 



1 



2(2tt) 



3 ' 



d 6 p 



0.078354 ±2 x 10" 
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= 0.235075 ±4 x 10" 6 . (53) 

According to Eq. (|5i~D , for fi < m^ n , (ijjip} = (ipip}^ ^ 0, i.e., chiral symmetry 

is spontaneously broken. For // > m^„, (ipip) = 0, i.e., chiral symmetry is restored. 
Therefore, there is a first order chiral phase transition and the critical value of fi is given 
by 

K = m &=p^> ( 54 ) 

The critical chemical potential \ic is equal to the dynamical quark mass at // = 0, which 
agrees with the result from an entirely different method [H| . (The authors argued this was 



a second order phase transition, in contrast we clearly observe a first order transition). Our 
result is consistent with other theoretical predictions He ~ because (see below) at 

fi = Q holds Mf « 3m£. 

2.3.4 Quark number density and susceptibility 

We can compute now the quark number density in the chiral limit m = 0, which yields 

-1 dE n (n\E x ^(x)^(x)\Q) 
n « = 2NJWs^~ W l = e(/i-„c), (55) 



which is consistent with the (3 = simulation results described in fll7| , and however, is 
different from the large \x behavior in the continuum (i.e. the Stefann-Boltzmann law 
n q oc /i 3 ). It remains to be seen whether higher order 1/g 2 calculations will improve this 
behavior. 

The quark number susceptibility, standing for the response of the quark number density 
to infinitesimal changes in /x, is 

dn q 
dfx 

2.3.5 Mass spectrum 



Xq = ^ = 5 (j* ~ f*c) • ( 56 ) 



Finally, let us look at some implications on the thermal mass spectrum of the pseudoscalar 
meson, vector meson and nucleon. The thermal mass is defined by Mfi = (h\H—fiN\h)—En. 
For the pseudoscalar meson, in the chiral limit m = 0, 

for ji < fi C , 



K = G, VI - MP = 0) = I for ; > - (57) 

Therefore, in the broken phase, the pseudoscalar is a Goldstone boson (M* oc \fm — > 0), 
and in the symmetric phase, it is no longer a Goldstone boson. For the vector meson, 

r(0) 



My = G lV /l-tanh 2 2^ = 0)= \ ^ < ^ (58) 

I 4m^ ; n for /i > /i C . 
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where My^ = 4y/d — l/(ag 2 CN) is the vector mass at fi = 0. Therefore, dM/d/i oc 
5(/j, — He) for the pseudoscalar and vector mesons. It is worth mentioning in Ref. fl8|], the 
authors found dM/dfi = outside the critical region. To see the critical behavior at zero 
temperature, one should be very close to \ic- This behavior is consistent with that of the 



quark number density discussed in Sect. [2.3.4| . To see whether the meson thermal masses 
depend on /i, higher order 1/g 2 corrections must be included. 
For the nucleon, we obtain the expected behavior 

M* = il4 0) -3^ (59) 
for /i < fic, where Mjy rs 3m^,. This leads to = at /i = fic- 



3 Outlook 

In this paper, we have developed a Hamiltonian approach to lattice QCD at finite density. 
It avoids the usual problem of either an incorrect continuum limit or a premature onset of 
the transition to nonzero quark density as fi is raised. The main result in the free case is 
given by Eq. ([34]), and those in the strong coupling regime are given by Eqs. (^0|)-([59|). 
We have seen that the approach works well in the free case and also in the strong coupling 
regime. We predict that at strong coupling, the chiral transition is of first order, and the 
critical chemical potential \ic ~ M^/3. 

Here we have only considered zero temperature. In the case of finite temperature, 
contributions from thermal excitations will make the calculations quite complicated. We 
plan to address this issue in a future paper. 

We are also aware that the strong coupling limit is not compatible with the continuum 
limit where a — > and 1/g 2 —>■ oo. For pure gauge theory, within a Hamiltonian approach, 
we can extend to the intermediate coupling and obtain meaningful results for the glue- 
balls [[HJ. For fermions, the calculation is far from trivial. Recently we proposed a Monte 
Carlo technique in the Hamiltonian formulation |2(| for the purpose to do nonperturbative 
numerical simulations, by combining the virtues of the Monte Carlo algorithm with impor- 
tance sampling and the Hamiltonian approach. We hope to apply it to QCD and with the 
aim to obtain useful information for RHIC physics. 
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